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Abstract 

Effective operators needed in single-beam two-photon transition calculations have been repre- 
sented with modified Goldstone diagrams similar to the type suggested by Duan and co-workers 
[J. Chem. Phys. 121, 5071 (2004) ]. The rules to evaluate these diagrams are different from those 
for effective Hamiltonian and one-photon transition operators. It is verified that the perturbation 
terms considered contain only connected diagrams and the evaluation rules are simplified and given 
explicitly. 
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I. INTRODUCTION 



Two-photon laser spectroscopy has been used as an important complementary technique 

in studying optical materials containing lanthanide ions. A recent short review can be found 

in Ref. Q. Two-photon excitation in a conventional microscope studying biological samples 

provides great advantages, including low photo-damage and depth discrimination, reduced 

photobleaching, good separation of the excitation light from the fluorescence emission. 

Two-photon lasers could give rise to entangled twin beams which offer potential for novel 
fl 

future applications. [3( A theoretical challenge in the field of two-photon spectroscopy is 



ilenge m 



prediction of two-photon transition rates. 

Our work focuses on the calculation of single-beam two-photon absorption rates, partic- 
ular for optical materials with lanthanides ions. In a former paper, [6j we have presented 
the effective operator method and given a perturbation expansion. In this paper we show 
that the perturbation expansion can be represented with a kind of diagrams, which can 
be calculated using one and two-particle matrix elements of interactions and single-particle 
energies. In Section II, we show how to derive the diagrams from the perturbation expan- 
sion; Section III presents those diagrams with only valence internal lines and unusual energy 
denominators; In Section IV we show the cancellation of disconnected diagrams. 



II. CONSTRUCTION OF DIAGRAMS FROM PERTURBATION EXPANSIONS 

The calculation of two-photon transition rates has been reduced to the calculation of an 
effective transition operator [6(. The square modulus of the matrix element of the effective 
operator between two eigenstates of the effective Hamiltonian (usually approximated with 
semi-classical phenomenological Hamiltonian) is the two-photon transition line strength that 
can be compared to measurements. In Ref. the effective operator is partitioned into 
two terms. The first term is calculated with "direct method" pi, H| , and the second term 
is calculated using many-body perturbation theory. It has been shown Ref. j^J that this 
method contains more general contributions than previous methods. Here we start from 
that expansion for the second term and consider the usual case where the model space Piq 
and Ppo are the same complete model space, denoted as Po- By "complete model space" 
we mean that the model space contains all the bases that can be formed by distributing 
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the valence electrons among the valence single particle states. 0, 0| Denoting the space 
orthogonal to Pq by Qo, the perturbative expansion is simplified into the following form 
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where T^ S i is the part of effective operator zeroth (i = 0) and first oder (i = 1) in V, k, 
ki and k2 sum over many-particle orthonormal eigenstates of H in Q Q and / sums over 
eigenstates of H in P . f and i are orthonormal eigenstates of H in P . 

The perturbation V contains both a one-body operator, denoted as V\ and a two-body 
operator, denoted as V 2 . Usually the operator O is one-body operator. Using algebraic 
representation and corresponding diagrammatic representation given in Ref|9j, Eq. 1 and 
Eq.4, O, V\ and Vi can be represented by 4, 4 and 9 diagrams respectively. Compositions of 
two O's and one V\ or V2, as given in Eq|2HEl with all possible contractions and all possible 



relative positions, resu 
diagrams given in Ref. 



t into hundreds of diagrams. This prompts us to adapt the modified 
to combine the set of diagrams with the same set of vertices, line 
directions and contractions but different relative position and internal orbit types into one 
diagram. 

A. Adapted connected diagrams with vertices V\ and two O 

There are altogether 3 ingoing and 3 outgoing lines, with at least 2 pairs of them con- 
tracted to form connected diagrams. Therefore the connected diagram can only be zero-body 
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(a constant) or one-body diagrams. A constant will not contribute to any transition, which 
is neglected throughout of this paper. 

There are three modified one-body diagrams as shown in Fig. HJ Similar to the the case 
for effective one-photon transition operators in Fig. 4 of Ref. |9|, each diagram actually 
represents contributions from several different terms in Eq|21-Eq|S| For example Fig. [TJa), 
with different internal line types, includes contributions from Eq|2J Eq|3]and EqlHl Note that 
the energy denominator rules will be slightly different from those for effective one-photon 
transition operators in RefQ 



B. Adapted connected diagram with vertices V2 and two O 

There are altogether 4 ingoing and 4 outgoing lines for for one Vi and two O. Therefore 
connected diagrams can be zero-body, one-body and two-body. 

There are also three modified one-body diagrams as shown in Fig. El , and five modified 
two-body diagrams as shown in Fig. El where two diagrams complex conjugate to FigOfb) 
and (c) are not given explicitly. 

Similar to connected diagrams in Fig^ Fig|2fa), includes contributions from Eq. 121 El 
13 These contributions form altogether 6 different normal diagrams. Note that the energy 
denominator rules for contributions from different terms are slightly different. 



C. Rules to evaluate these diagrams 

Remembering that each modified diagram actually represents a set of diagrams with 
different ordering of vertices and different types of internal lines, it is not difficult to obtain 
the analytic expressions for these diagrams. For each diagram with fixed vertex ordering, the 
expression can be calculated by the same rules 1)— 7) given in Sec. II, Ref. [n|, except that 
the rule to calculate energy denominators in rule 6) should be modified as follows: Denoting 
the net outflow energy of a loop enclosing vertices 1, 2, • ■ ■ , % as 

E noe (l, 2, ■ ■ ■ , i) = Koe(l) + £noc(2) + • • • + £U(i), (9) 
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and supposing the number of O in the loop is io {io = 0, 1, 2) and the number of vertices 
for the diagram is n, the energy denominator can be written as 

D = IS (y£noe(l, 2, ■ • • , n) - E noe (l, 2, ■ ■ ■ , i)) . (10) 

III. DIAGRAMS WITH VERTEX OPERATORS ACTING ONLY ON MODEL 
SPACES P 

Eq. 01 and |H1 give diagrams with some vertices, whose bra states and ket states are 
in model spaces Pq. As an example, Eq. HJhave the same operator sequence OVO as that 
of Eq. |U1 but different denominators. They will give diagrams with different evaluation 
rules, similar to but not exactly the same as those terms in effective Hamiltonian that give 



folded diagrams. 



n 



12, 



131 ] Instead of introducing new folded diagrams, we combine these 
contributions to the diagrams given above. The rules to evaluate these contributions differ 
from standard rules in an overall factor —1/2 for diagrams from EqEJ and EqlHl and in 
energy denominators by using rules given in Eq0J EqEJ Eq0 and EqlHl 



IV. CANCELLATION OF DISCONNECTED DIAGRAMS 

For modified disconnected diagrams with one V vertex and two O vertices, there must 
be a connected part with only one vertex, which is either V oi O. 

Case A: V is disconnected from other parts. In this case both the ingoing and the 
outgoing lines of V are valence lines. Such a diagram includes contributions from EqEHHl 
with the same matrix elements but different coefficients due to different energy denominators. 
Suppose the net-outgoing-energies from the first O and the second O vertices are AE\ and 
AE2 respectively, and that from the V is A £3. The three different coefficients for the three 
contributions from EqEHU respectively, follows 



1 



■ [AE 1+ AE 3 +Ai; 2 )+0 _ + AE ^ (AE 1+ AE 3 +AE 2 )+0 _ ^ 



(11) 



1 1 

+ (-5) [( (AE 1+ AE 2+ AE 3 )+0 _ ^ ~ AE ^AE 1+ AE 2 +AE 3 )+AE 3 _ + ( 13 ) 
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(14) 



Case B: one O is disconnected from other parts. Such a diagram includes contributions from 
Eq. I2HI3 It can be verified in the same way as above that all these contributions cancel. 

V. CONCLUSION 

Diagrams for two-photon transitions rates have been constructed and have been shown 
to be connected. Simplified one-photon diagrams given in J. Chem. Phys. 121, 5071 (2004) 
have been adapted to two-photon cases. The algebraic expressions suitable for a general 
two-photon calculation can be given from these diagrams. 
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FIG. 1: Modified one-body diagrams first order in V\ for one-beam two-photon transition moment. 
Note that: 1) the internal lines can be both particle lines and core lines, and the relative altitudes 
of two O's and V\ change accordingly; 2) contributions of (b) and (c) to the effective operator are 
complex conjugation of each other. 
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FIG. 2: Modified one-body diagrams first order in V2 for one-beam two-photon transition moment. 
Note that: 1) the internal lines can be both particle lines and core lines, and the relative altitudes 
of two O's and V\ change accordingly; 2) there are natural constrains over several internal lines 
forming a close loop, since particle lines go upward and core lines go downward, for example, in (b), 
{j, k} must be one particle line and one hole line; 3) contributions of (b) and (c) to the effective 
operator are complex conjugation of each other. 
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FIG. 3: Modified two-body diagrams first order in V2 for one-beam two- photon transition moment. 
Note that: 1) the internal lines can be both particle lines and core lines, and the relative altitudes 
of two O's and V\ change accordingly; 2) two diagrams with contributions complex conjugate to 
these of (b) and ( c) are not shown here. 
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